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^"*^ ' A shear-improved Smagorinsky model is i ntroduced based on recent results concerning 

shear effects in wall-bounded turbulence bv lToschi et all (|2000|) . The Smagorinsky eddy- 
viscosity is modified as vt — (C s A) 2 (|5| — | (5) |): the magnitude of the mean shear 1(5)1 is 
subtracted to the magnitude of the instantaneous resolved strain-rate tensor |5|; here C<? 
is the standard Smagorinsky constant and A denotes the grid spacing. This subgrid-scale 
model is tested in large-eddy simulations of plane-channel flows at Reynolds numbers 
Re T = 395 and Re T — 590. First comparisons with the dynamic Smagorinsky model 
and direct numerical simulations, including mean velocity, turbulent kinetic energy and 
Reynolds stress profiles, are shown to be extremely satisfactory. The proposed model, in 
\f~^ addition of being physically sound, has a low computational cost and possesses a high 

potentiality of generalization to more complex non-homogeneous turbulent flows. 
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1. Introduction 



The prohibitive cost of direct numerical simulations (Dns) of turbulent engineering 
flows motivate the elaboration of simplified models, requiring less com putation but still 
j> \ relevant (to some deg r ee) for reproducing the lar ge-scale dynamics ()Deardorfll Il97flt 

■ iLesieur fc Metaislll996t IPiomellilll999t ISagautll200ll^ . In this context, the modeling near 

r> \ a solid boundary is of particular interest. A boundary affects the kinetics of the flow 

through different mechanisms; the most prominent is that related to the mean shear, 
which is extremal in the proximity of the boundary and responsible for the production 
of streamwise vortices and streaky structures. These flui d structures eventua lly detach 
from the boundary and sustain turbulence in the bulk l)Perot fc Moinlll995h . Thus, it 
is thought that understanding how the mean shear impacts on fluid motions is a key 
to improving the capabilities of numerical models of turbulent flows. Turbulence near 
a boundary is not homogeneous nor isotropic and therefore, the customary theor etical 
backg round for modeling the unresolved small-scale dynamics should be surpassed (Pope 
12000^1 . In the present work, we present a new subgrid-scale model of turbulence which 
may be viewed as an improvement over the popular Smagorinsky model in the presence 
of a mean shear. Our model originates from recent theoret ical findings concerning shear 
effects in wall-bounded turbulence, sec iToschi et al\ l|2000j) . We propose to subtract the 
magnitude of the mean shear from the resolved strain-rate tensor in the definition of 
the eddy- viscosity. As this will be discussed later, this simple improvement accounts for 
the non-isotropic nature of the flow near the boundary and at the same time, allows 
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us to recover the standard Smagorinsky model in regions of (locally) homogeneous and 
isotropic turbulence. The general framework of the so-called large-eddy simulation (Les) 
of turbulent flows is now briefly recalled. 

Roughly speaking, large-scale motions transport most of the kinetic energy of the flow. 
Their strength make them the most efficient carriers of conserved quantities (momentum, 
heat, mass, etc.). On the contrary, small-scale motions are primarily responsible for the 
dissipation while they are weaker and contribute little to transport. From mechanical 
aspects, the large-scale (energy-carrying) dynamics are thus of particular importance and 
the costly computation of small-scale dynamics should be avoided. Furthermore, while 
large-scale motions are strongly dependent on the external flow conditions, small-scale 
motions are expected to behave more universally. Hence, there is a hope that numerical 
modeling can be feasible and/or require few adjustments when applied to various flows. 

In Les, only the large-scale components of flow variables are explicitly integrated in 
time; interactions with the unresolved small-scale components are modeled. A spatial 
filtering is conceptually introduced as 0(x, t) = J </>(x', £)Ga(x — x')cZx', where the filter 
width A fixes t he size of the smallest scales of variation retained in the flow variable 0(x, t) 
<)Leonardll974l) . In practical applications, A is chosen much larger than the spatial cutoff 
scale of 4>(x,t), i.e. the dissipative scale of turbulence, so that </>(x,t) may be considered 
as the large-scale component of 0(x, t). Applying the previous filtering procedure to the 
Navie r-Stokes equations (and neglecting here non-commutation errors l|Ghosal fc MoirJ 
Il995|) for the sake of simplicity) yields 

dui | _ dig | dTjj _ ^P_ +v d 2 "» with duj_ _ Q - ^ 

dt 3 dxj dxj dxi dxkdxk dxi 

where Ui (x, t) and J>(x, t) represent the large-scale velocity and pressure, and v is the kine- 
matic viscosity of the fluid. The equations 1|1.1[) arc amenable to numerical discretisation 
with a grid spacing comparable to A, since u(x, t) is expected to vary smoothly over the 
distance A. Tij(x,t) = Ui(x,t)uj(x,t) —Ui(x.,t)uj{x.,i) is named the subgrid-scale (Sgs) 
stress tensor and encompasses all interactions between the grid-scale component and the 
unresolved subgrid-scale component of u(x, i). In Les, Tij(x, t) needs to be e xpressed 
in ter ms of the grid-scale velocity field u(x, t) only, which is the hard problem ijLesieurl 
Il997j) . Eddy- viscosity models parameterize the SGS stress tensor as 

Uj - -SijTkk = -2v T Sij where (x,t) = - ^^-(x,t) + -^-(x,t)j (1.2) 

where i*r(x, t) is the scalar eddy-viscosity and Sij is the resolved rate-of-strain tensor. 
This empirical modelization is rooted in the idea that SGS motions are primarily re- 
sponsible for a diffusive transport of momentum from the rapid to the slow grid-scale 
flow-regions. The theoretical basis for the introduction of an eddy- vis cosity is rather ins e- 
cure, however, it appears to be workable in practice (as advocated bv lKraichnanl l)l976|) '). 
The eddy-viscosity vt is then primarily designed to ensure the correct mean drain of 
kinetic energy from the grid-scale flow to the SGS motions: —(Tij Sij) from the equations 
Ql.lfl . Another important feature is that ^r(x, t) should vanish in lam inar flow-regions, 
e.g. in the viscous sub-layer near the boundary l|Moin fc Kimlll982l) . Accordingly, our 
main concern was to determine an eddy- viscosity vt (x, t) that would take into account 
shear effects in the transfer of kinetic energy to the SGS motions (without any sort of 
dynamical adjustment) and naturally decrease to zero at the boundary (without using 
any ad-hoc damping function). 

In Sectional the advantages and drawbacks of the Smagorinsky model, which is used as 
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our baseline model, are recalled. The necessity for a new model that can achieve a satis- 
factory compromise between accuracy and manageability is identified. A shear-improved 
Smagorinsky model is introduced within the classical picture of shear turbulence. In Sec- 
tion |31 we present results from Les of turbulent plane-channel flows at Re T = 395 and 
Re T = 590. Comparisons are carried out with the dynamic Smagorinsky model and Dns. 
Discussion and perspectives follow in Section 0] 



2. A shear-improved Smagorinsky model 

2.1. Our baseline model: The Smagorinsky model 

The Smag orinsky model i s cert ai nly the simp lest and most commonly used eddy- viscosity 
model ( see ISmaeorinskvl (|l963j) . |popel l)2000(l for a comprehensive description). The pre- 
scription for vt writes i*r(x, i) = (C s A) 2 |S'(x, t)\, where \S\ = {2SijSij) 1 / 2 represents the 
magnitude of the resolved rate-of-strain and C s is a non-dimensional coefficient (called the 
Smagorinsky constant). The major merits of the Smagorinsky model are its manageabil- 
ity, its computational stability and the simplicity of its formulation (involving only one 
adjusted parameter). A ll this makes it a very valuable tool for engineering applications 
l|Rogallo fc Moinlll984|) . However, while this model is found to give acceptable resu lts in 
the L es of homogeneous and isotropic turbulence (with C s « 0.17 according to l)Lillvl 
it is found to be too dissipative with respect to the resolved motion s in near- wall 
turbu lence, due to an excessive eddy- viscosity arising from the mean shear 1 Moin fc Kiml 



Il982l) . The eddy- viscosity predicted by Smagorinsky is non-zero in laminar flow- regions; 
the model introduces spurious dissipation which has the effect of da mping the growth 
of sm all perturbations and thus restrain the transition to turbulence l|Piomelli fc Zand" 
Il99lh . 

To alleviate these deficiencies in the case of wall-bounded flows, the Smagorinsky con- 
stant C s is often multiplied by a damping factor de pending on the w all-normal distance, 
the van Driest function being the prime example Ijvan DriestHl95fift . Although the van 
Driest damping function is commonly employed, its theoretical basis has never been 
adequately addressed, thereby leaving this function rather arbitrary. Moreover, if the de- 
termination of the distance is straightforward in the case of a plane boundary, it becomes 
more ambiguous near a curved boundary or a sharp corner. A more general approach, 
free from the use of the wall-normal distance, is therefore desirable. The dynamic SGS 
model intends to evalu ate the Smagorinsky c onstant (from the resolved motions) as the 
calculation progresses I Germano et oiJfl99ll) and thus avoids the need to specify a pri- 



ori, or tune, its value. In brief, the adjustment of C s relies on the Germano identity 
l)Germa,nolll99l and as sumes the scale similari ty of resolved velocity fluctuations at 
scales comparable to A (|Meneveau k, Kat3l2000|) . This methodology yields a coefficient 
C s that varies with positio n and time an d vanishes in the vicinity of a solid boundary 
with the correct behavior l|Piomellilll9'9l . It is beyond dispute that the dynamic pro- 
cedure greatly improves the capacity of the original Smagorinsky model, however, this 
progress is also accompanied with a certain deterioration of the numerical stability and 
an increase of the computational cost. In this situation, we believe that it is legitimate to 
seek for a SGS model that would achieve a better compromise between accuracy and man- 
ageability: ideally, as simple as the original Smagorinsky model and as accurate as the 
dynamic Smagorinsky model. The present SGS model that we shall call shear-improved 
Smagorinsky model has been elaborated along this line of idea. 
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2.2. Our proposal: A shear-improved eddy-viscosity 

Our SGS model differs from the standard Smag orinsky model in h ow the e ddy- viscosity is 
define d. It borrows ideas originally advanced bv lSchumannl jl975(l (see also lSullivan et al\ 
<jl994^) 1 by using a two-part eddy- viscosity accounting for the interplay between the non- 
linear energy cascade present in isotropic turbulence, and mean shear effects associated 
with anisotropy. However, it differs from Schumann's proposal in its formulation, in 
particular for an extra simplicity. Our proposal for the SGS viscosity writes 

ut = (C s A) 2 .(\S\-\(S}\) (2.1) 

where the brackets ( } would a priori denote an ensemble average, in practice, space 
average over homogeneous directions and/or time average will be used (this specific issue 
will be mentioned again in Section^}. From the definition 12.111 the mean drain of kinetic 
energy from the grid-scale to the subgrid-scale motions is given by 

F sgs = -(nM = (C s Af ((\S\ 3 ) \(S)\{\S\ 2 )) ■ (2.2) 

Straightforwardly, our eddy-viscosity vanishes if the resolved turbulence disappears, i.e. 
if S = (S). We shall argue that it is also consistent, to some extent, with the SGS energy 
budget of sh e ar tur bulence. The theoretical basis of our model was first put forward by 
iToschi et all l)2000|) . on account of previous numerical and experimental studies on wal l- 
bounded turbulence ^Toschi et ollll999l iBenzi et al.lll999l iRuiz-Chavarria et all\2()0(^ . 
For the sake of simplicity, we shall here recast t he key arguments in th e ideal case of a 
statistically stationary homogeneous shear flow I Monin fc Yaglorr]ll975l) . 



In an homogeneous shear flow, the velocity field u(x, t) may be decomposed into 
iti(x, t) = ?4(x, t) + (dUifdxj)xj , where Ui(x) and m^(x, t) denote the mean and fluc- 
tuating part of the velocity, respectively. Starting from the exact dynamical equations 
for the two-point correlation function R(r) = (u£(x, i)uj(x + r, t)) and by integrating this 
equation over a sphere B r of radius r centered at x one ca n establish an energy budget 
(at scale r) which takes the form fsee lCasciola et all (|2003^ l : iDanaila et 

5r(r) + 5 3 pr (r) = -ier + 2u± j (|<5u'(x, r, t)| 2 ) ds) . (2.3) 

The two contributions S^(r) and S^ r (r) arise from the non-linear term of the Navier- 
Stokes equations. In the right-hand side of H2.3|) e denotes the mean energy dissipa- 
tion rate and the second term encompasses finite Reynolds effects (at scale r) with 
5u'(x, r, t) = u'(x + r,t) — u'(x, t). The energy budget 12.3fl is the general ization of th e 
Karman-Howarth equation in the context of homogeneous shear turbulence l)Hinzelll976^ . 
In the context of Les, it may also be interpreted as the proper mean SGS energy budget 
with respect to the grid scale r. More explicitly, 



5 3 r W = ^f dB (<l<5u'(x,M)| 2 <K(x,r,*)} + gr.-flJu'^M)! 2 )) dS, 



(2.4) 



represents the transfer of kinetic energy from grid-scale motions (at scales larger than r) 
to subgrid-scale motions. Also, it indicates that this transfer results from both the grid- 
scale turbulent fluctuations and the mean shear. These two effects correspond respectively 
to the non-linear triple-correlation term and the rapid (or linear) t erm entering in the 
spectral decomposition of the energy transfer, as first evidenced bv ICraval l)l95S|) . The 
second term in the left-hand side of (|2.3|) represents a production of SGS kinetic energy 
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induced by the mean shear, and is expressed as 

SfM = J B 2^(<y«i(x,r,*)^(x,r,t)>dV. (2.5) 

In summary, the energy budget (|2.3|) means that in the steady-state, the SGS dynamics 
are sustained against molecular dissipation (represented by the right-hand side of (12.311 ) 
by the transfer of energy originated from grid-scale motions (S\ T (r)) and the production 
of energy induced directly by the mean shear (S^ r (r)). 

As mentioned in the introduction, the estimate of the SGS energy flux, which should 
be identified with Srf(A)/A in the previous developments, is of prime importance in the 
modeling of the eddy- viscosity. Roughly speaking, one may consider that turbulent grid- 
scale velocity differences 5u'(A) typically behave as the fluctuating part of the resolved 
rate-of-strain [S 7 ] multiplied by A, i.e. Su'(A) ps IS 7 ] A. In the same way, SU(A) w \ (S)\A. 
Equation (|2.4(l then suggests that the mean SGS energy flux should involve two separate 
contributions of order A 2 jJS T \ 3 ) and A 2 |(5}|(|5 7 | 2 ), respectively. 

In flow-regions where \S'\ 3> \{S)\, the SGS energy flux reduces to the contribution of 
order A 2 (|S''| 3 ). In these regions, the mean shear is too weak to perturb the grid-scale 
dynamics; eddies of size comparable to the grid-scale A adjust dynamically via non- 
linear interactions to transfer energy to SGS motions. This is the standard mechanism 
behind homogeneous and isotropic turbulence llFrischlll995|) . From our expression (|2.2() . 
we remark that the Smagorinsky estimate F sgs ~ (C S A) 2 (\S'\ 3 ) is consistently recovered 
in that case. In regions where \(S)\ ^> \S'\, the behavior of the flow is clearly different: 
eddies of size comparable to the grid -scale have no t ime to adjust dynamically and are 
rapidly distorted by the mean shear ijLiu et aZ.lll999|) . In these regions, the SGS energy 
flux is driven by the mean shear and therefore is dominated by the contribution of 
order A 2 |(5)|(|5 7 | 2 ). From expression lj2~2~)> and assuming that (\S\ 3 ) « (\S\ 2 ) 3/2 , one 
obtains F sgs ~ 1/2 (C S A) 2 \(S)\ (IS'I 2 ) in agreement with the previous reasoning. We 
may thus conclude that proposal (|2.1J) for the eddy-viscosity is consistent with the SGS 
energy budget of (locally homogeneous) shear turbulence in the two limiting situations 
IS"' | 3> \(S) \ and \ {S)\ \S'\. The first results concerning plane-channel flows (see Section 
[3J indicate that the proposed model abridges between these two situations without the 
need for any additional adjustment. 



3. Les of turbulent plane-channel flows 

Over the l ast twenty years Le s of w all-bounded flows have received a considerable 
attention (see IPiomelli fe Balarasl l|2002() for a recent review), with the turbulent plane- 
channel flow (Kim"" eFoITijl987(l for an original work) being the prototypical case. This 
flow allows for the investigation of shear effects in a simple geometry and has therefore 
provided a useful test bed of our eddy- viscosity model. Furthermore, we have been able to 
confront our results, including mean velocity, turbulence intensities and Reynolds stress 
profiles, with the well-estab lished literature pre sent on that case, e.g. the comprehensive 
Dns d atabase obtained by iMoser et al\ ljl999|) or, more recently, by Irlovas fe Jimene3 

3.1. Numerical simulations 

We performed two Les at Re T = 395 and Re T — 590, where Re T is the Reynolds 
number based on the friction velocity u T : Re T = u T H/v (H is t he half width of th e 
channel) . These two cases correspond to the Dns conducted by IMoser et all ()l999l) . 
Periodic boundary conditions were imposed in the streamwise and spanwise directions; 
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asymptotic mean Re =590 (96x97x96) 
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Figure 1. The time development of the Reynolds number Re T (t) based on the plane-averaged 

friction velocity u T (t). 

no-slip conditions at the wall. The flow was simulated by integrating the filtered Navier- 
Stokcs equations with the prescribed eddy- viscosity i|2.1[) . The right-hand side of 

the equations i|l.lf) was supplemented by an (external) pressure-gradient Ap ext /47ri? in 
order to drive the flow in the streamwise direction. The pressure difference Ap axt (t) was 
adjusted dynamically to keep a constant flow rate through the channel. The integration 
in time relied on a Fourier-Chebyshev pseudo-spectral solver (de-aliased by using the 
3/2 rule) based on a third-o rder Runge-Kut ta scheme. More details about the numerical 
simulations can be found in lXu et all ([l996'l. The Smagorinsky constant was fixed to its 
standard value C s = 0.16 (in the case of homogeneous and isotropic turbulence) and the 
scale A was estimated as (AxAyAz) 1 / 3 , where Ax, Ay and Az denote the local grid 
sizes in each direction. In the following, the grid-scale velocity components arc U + u' , v' 
and w' along the streamwise, wall-normal and spanwisc directions, respectively. 

3.2. Numerical results 

At initial time, velocity distributions were designed to satisfy the Poiseuille profile plus 
a small random perturbation. The time development of Re T (t) = u T (t)H / 'v , where u T (t) 
expresses as the square root of the plane-averaged wall shear stress, is plotted in Figure 
^ for the case Re T = 590. A transition (or drag crisis) to the appropriate turbulent 
regime occurs naturally as the integration (involving the eddy-viscosity (|2.1[) ') progresses. 
This feature constitutes a first improvement over the Smagorinsky model, for which such 
transition is not captured. Once a developed turbulent state was achieved, statistics were 
accumulated. 

The mean velocity U + (y + ) is displayed as a function of the wall- normal distance 
y + in Figure From now on, the average is meant in time and over horizontal x-z 
planes (homogeneous directi ons). Mean velocity profiles agree very well with the Dns 
data of reference obtained bv lMoser et al\ l)l999j) for both considered Reynolds numbers. 
The c omparison with the dynamic Smagorinsky model is also very satisfactory ijPiomellil 
1993|) . Turbulence intensity profiles (normalized by the squared friction velocity) are 
displayed in Figure Also here, the profiles compare very well with the Dns data 
and les results based on the dynamic Smagorinsky model: notice how peak values are 
suitably captured. As for the dynamic Smagorinsky model, a slight drop is observed 
for the spanwise and wall-normal components, in the log-layer; on the contrary, the 
streamwise component fits remarkably the Dns data. In the proximity of the wall, the 
correct behavior is obtained. 





1Cf 1 10° 10' 10 2 10"' 10° 10' 10 2 



y + y + 

Figure 2. Left: (•) mean-velocity profile (in wall units) at Re T = 395. The computational 
domain (in outer units) is 4tvH x 2H x 2nH w i th 64 x 65 x 64 grid points. In comparison 
with (— ) the Dns data obtained bv lMoser et al\ (I1999T) in the domain 2-kH x 2H x irH with 
256 x 193 x 192 grid points, and (A) a computation of the dynamic Smagorinsky model carried 
out by Piomelli (private communication) in the domain 5nH/2 x 2H x nH/2 with 48 x 49 x 48 
grid points (using a pseudo-spectral solver). Right: (•) mean- velocity profile at Re T = 590 with 
96 x 97 x 96 grid points. In comparisons with (— ) the Dns data with 384 x 257 x 384 grid points. 

Re =395 (64x65x64) Re =590 (96x97x96) 




y + y + 



Figure 3. Left: Turbulent intensity profiles at Re T — 395 in comparison with Dns data ( — ) and 
Les data obtained with the dynamic Smagorinsky SGS model (open symbols). Right: Turbulent 
intensity profiles at Re T = 590. The inset focuses on the near-wall behavior in comparison with 
Dns data. 



These positive results make us presume that our eddy-viscosity has captured the es- 
sential of shear effects in wall-bounded turbulence. Now, in order to better check this 
hypothesis, the ratio y ( | S' | 2 ) / is shown as a function of the wall-normal distance 
in Figure For y + < 25, one can see that the mean shear dominates over the fluctuat- 
ing part of the rate-of-strain, indicating the predominance of the mean-shear component 
of the eddy-viscosity in that region. Note that the transition distance y + ~ 25 is fully 
consistent with the empir ical charact eristic distance A + = 25 commonly used in the van 
Driest damping function l|Ponell200nl) . In the log-layer, ^(\S'\ 2 ) / \(S)\ increases slowly 
with y + (the standard descriptio n of the log-layer p redicts a linear increase, according to 
oc 1/y and (IS 1 '! 2 } oc u 2 /A 2 (|Schlichtindll95. t ^ ) and eventually diverges around the 
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Re T =395 (64x65x64) 
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Figure 4. The ratio \J (IS*'! 2 ) / \(S)\ is displayed as a function of the wall-normal distance y + 
in the Les at Re T = 395. For y + < 25, the mean shear dominates over the fluctuating part 
of rate-of-strain. In that region, our eddy-viscosity differs from the original Smagorinsky model 
candidate. 



Re_=395 (64x65x64) Fte T =395 (64x65x64) 




1(T 1 10° 10 1 10 z 50 100 150 200 



Figure 5. Left: The viscous and Reynolds stresses (computed from the resolved velocity) at 
Re T — 395. The two contributions are equal at y + ~ 12, in full agreement with Dns results. 
Right: The peak production of turbulent energy also occurs at y + ~ 12, as expected. 



ccntcrlinc of the channel. Thus, wc may claim that our model suitably abridges between 
the situation where the mean shear prevails (close to the boundaries) and the situation 
where the fluctuating part of the rate-of-strain dominates (in the bulk of the channel) . 

The viscous, dll + /dy + , and the Reynolds. (u + v + ), contributions to the resolved stress 
are shown i n Figure [S] They equal for y + « 12 in perfect agreement with Dns results 
llPopel200(t . The peak value of the turbulent energy production, (u + v + )dU + /dy + , occurs 
at the same distance y + w 12, as expected theoretically from the Navier-Stokes equations. 
The grid-scale Reynolds stress profile is displayed in Figure [H] The agreement with the 
Dns and the Les based on the dynamic Smagorinsky model is fair and also the behavior 
close to the wall is very satisfactory. 
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4. Discussions and perspectives 

Concluding, we presented a very simple SGS model consisting in a physically sound im- 
provement over the well-established Smagorinsky model. An explicit connection with the 
scale-by-scale energy budget of homogeneous shear turbulence has also been identified. 
Our first results for turbulent plane-channel flows indicate that the proposed model pos- 
sesses a very good predictive capability (essentially equivalent to the dynamic Smagorin- 
sky model) with a computational cost and a manageability comparable to the original 
Smagorinsky model. 

The generalization to more complex non-homogeneous flows is a priori straightforward 
since no geometrical argument enters in the definition of the eddy-viscosity. However, 
an appropriate average must be specified in the absence of homogeneity directions. As 
a natural candidate, we suggest an average in time to evaluate the mean components of 
the rate-of-strain tensor (assuming that the flow is statistically stationary). In the case 
of non-stationary flows, an ensemble average (over several realizations of the flow) may 
be envisaged. These points are the subject of current investigations in progress. 

The authors are grateful to Ugo Piomelli for providing numerical data; in particular 
F.T. wishes to thank him for the many useful discussions. Fang Le helped us to perform 
the Les at the Tsinghua University in China. This work was jointly supported by the 
Ecole normale superieure de Lyon and the Ecole centralc dc Lyon under a Bonus-Qualite- 
Recherche grant. Finally, F.T. acknowledges E.L. for his kind hospitality at the Ens- 
Lyon. 
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